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Abstract 

We present a common chiral power-counting scheme for vector, axial-vector, scalar, and pseudoscalar WIMP-nucleon interactions, 
and derive all one- and two-body currents up to third order in the chiral expansion. Matching our amplitudes to non-relativistic 
effective field theory, we find that chiral symmetry predicts a hierarchy amongst the non-relativistic operators. Moreover, we 
identify interaction channels where two-body currents that previously have not been accounted for become relevant. 
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1. Introduction 

Elucidating the nature of dark matter is one of the most 
pressing challenges in contemporary particle physics and astro¬ 
physics. Still, one of the dominant paradigms rests on a weakly- 
interacting massive particle (WIMP), such as the neutralino 
in supersymmetric extensions of the standard model (SM). A 
WIMP can be searched for at colliders, in annihilation signals, 
or in direct-detection experiments, where the recoil energy de¬ 
posited when the WIMP scatters off nuclei is measured. Recent 
years have witnessed an impressive increase in sensitivity, e.g., 
from XENON 100 fl], LUX H, and SuperCMDS H, which 
will further improve dramatically with the advent of ton-scale 
detectors, XENON IT |01 and LZ ||3l. In the absence of a sig¬ 
nal, direct-detection experiments provide more and more strin¬ 
gent constraints on the parameter space of WIMP candidates. 
To derive these constraints and to interpret a future signal, it 
is mandatory that the nucleon matrix elements and the nuclear 
structure factors, which are required when transitioning from 
the SM to the nucleon to the nucleus level, be calculated sys¬ 
tematically and incorporate what we know about QCD. 

Effects at the level of the nucleus can be described by an 
effective field theory (EET) whose degrees of freedom are non- 
relativistic (NR) nucleon and WIMP fields |0|3]- This NREET 
has been recently used in an analysis of direct-detection exper¬ 
iments iQ]. In this approach, scales related to the spontaneous 
breaking of chiral symmetry of QCD are integrated out, with 
the corresponding effects subsumed into the coefficients of the 
EET. In the context of nuclear forces, such an EET is called 
pionless EET. To derive limits on the WIMP parameter space, 
information from QCD has then to be included in the analysis 
in a second step. 

Alternatively, one can start directly from chiral EET (ChEET) 
to incorporate the QCD constraints from chiral symmetry |9l- 
which makes predictions for the hierarchy among one- 
and two-body currents. Based on ChEET, scalar and axial- 
vector two-body currents were recently considered in 


and ifTllfl^ . respectively. Moreover, lattice QCD can be used 
to constrain the couplings of two-body currents iB. 

The goal of this Letter is to combine vector, axial-vector, 
scalar, and pseudoscalar interactions in a common chiral power 
counting, collect all relevant one- and two-body matrix ele¬ 
ments, and match the result onto NREET. This combines our 
knowledge of QCD at low energies; the one-body matrix ele¬ 
ments correspond to the standard decomposition into form fac¬ 
tors, while the two-body scalar lH, [IS, vector lITsl - liol] . and 
axial-vector EEl] currents have been calculated as well, 
the vector current even at one-loop order. Here, we combine 
these results for their application in direct detection, extending 
the axial-vector two-body currents to finite momentum trans¬ 
fer and generalizing to the three-flavor case where appropriate. 
By matching to the NREET, we find that the chiral symmetry 
of QCD predicts a hierarchy among the different operators and 
that two-body currents can be as important as one-body currents 
in some channels. 


2. Effective Lagrangian and kinematics 

We start from the following dimension-6 and -7 effective La¬ 
grangian for the interaction of the WIMP x, assumed to be a 
SM singlet, with the SM fields 1221 


A = ^ X xijsX mqqq 


+ Cfxxniqqijsq + CfxiysXmqqiysq] 

+ ^ X fVxyA:qy^iq + cfxyysxqy^q 
+ c^'^xY'x qy/^ysq + cfxy^'ysxqyiiysq] 

+ ^ X + Cfxo^''iy5Xqo-i,yq\ 







+ j^[clxxasG%G>: + Clxiy5XasG%G>: 

+ ClxxasG%&: + ClxiysX<ysG%&:\ (1) 

where the Wilson coefficients C, parameterize the effect of new 
physics associated with the scale A (organizing the interactions 
in this way assumes A to be much larger than the typical QCD 
scale of 1 GeV). To render the scalar and pseudoscalar matrix 
elements renormalization-scale invariant we included explicitly 
the quark masses niq in the definition of tfie respective opera¬ 
tors. We further assumed to be a Dirac fermion (in the Majo- 
rana case, — Cj/ — 0), and defined the dual field 

strength tensor as 




( 2 ) 


with sign convention = -i-l. Compared to the operator 
basis used in 02311 we do not include the dimension-8 opera¬ 
tors related to the traceless part of the QCD energy-momentum 
tensor. As shown in ll23n . these operators become relevant for 
heavy WIMPs and contribute to spin-independent interactions, 
decreasing significantly the single-nucleon contribution. Fi¬ 
nally, we will ignore the tensor operators in O and concentrate 
on the chiral predictions for the V,A,S,P channels. 

The kinematics for the WIMP-nucleon scattering process are 
taken as 

N{p)+x(k)^N{p’)+x(k'), (3) 


the momentum transfer is defined as 


q - k' - k - p- p', - t, (4) 


and the pion, q, nucleon, nucleus, and WIMP masses will be de¬ 
noted by M^, M,j, niff, niA, and m^, respectively (Dirac spinors 
are normalized to 1). We will also need 


P^p + p', K = k + k'. (5) 


The cross section differential with respect to momentum 
transfer for the elastic WIMP-nucleus scattering process in the 
laboratory frame can be expressed as 


dcr 

dq2 


1 

87tv\2J+ 1) 


2|MNRl'+0(q‘’), 

spins 


(6) 


with nucleus spin J, WIMP velocity v, and NR amplitude AInr 
defined as 

M - + 0{q^), (7) 

where M is the relativistic scattering amplitude. In the Majo- 
rana case, ® receives an additional factor of 4. 


3. Chiral power counting 

We use the standard chiral power counting 

d = 0(p), nig = O(p^), Ufj, Vfj = 0(p), (8) 

with axial-vector and vector sources and v^. The velocity 
distribution in dark matter halo models indeed suggests to count 


the momentum transfer q < M„ as 0(p) II^l. In the b ary on 
sector we depart from the standard counting in chiral pertur¬ 
bation theory (ChPT) and adopt the more conventional ChEFT 
assumption (see, e.g., |2^28]) for the scaling of relativistic cor¬ 
rections 


niN 


(9) 


This counting is appropriate for a break-down scale around 
500 MeV. As far as the WIMP is concerned, a chiral count¬ 
ing is only required for the NR expansion of the spinors. We 
assume the same counting as in the nucleon case, but display 
the corresponding additional powers explicitly. If > niff, 
the suppression will be more pronounced, for M„ ^ ni^ ^ niff 
the counting should be adapted, and for even smaller ni^ the 
naive counting breaks down. 

For most of the channels it suffices to consider tfie leading- 
order Lagrangian to determine at wfiich chiral order a given 
contribution starts. For the one-body matrix elements higher 
orders are subsumed into the nucleon form factors, which are 
obtained by their chiral expansion or could be taken from phe¬ 
nomenology. In this work, we consider all contributions up 
to 0{p^). Since the leading two-body terms start at O(p^), 
this leaves the possibility that the next-to-leading-order (NLO) 
pion-nucleon Lagrangian involving the low-energy constants 
Ci lE^ could be required, and this is indeed the case for the spa¬ 
tial component of the axial-vector current 111 lL|l2|] (indicated by 
“2b NLO” in Table [T]i. In the same channel, NN contact terms 
di lE^ enter. We define both c,- and di in the conventions of 10 
(with dimensionless cg and c?). 

As a preview of our results, the leading chiral orders of one- 
and two-body currents for time and space components of the 
axial-vector and vector currents, as well as for the scalar and 
pseudoscalar operators, are fisted in Table [T] The suppression 
by two powers (“-I-2”) originating from the WIMP spinors is 
displayed separately. In the following sections, we give results 
for all one- and two-body currents involved in Table [T] 


4. Nuclear matrix elements 

4.1. Scalar 

At zero momentum transfer the scalar couplings of the heavy 
quarks Q - c,b,t can be determined from the trace anomaly of 
the QCD energy-momentum tensor 0 




^QCD 


{N^N) = lUff, 


Pqcd _ ^ _ 2Nf Us 


11 


+ <9( q ' j ). 


2gs \ 3 I8n 

For Nf -2 active flavors, one obtains 

{N\mQQQ\N) = -^<A|G;,Gf |A) = 

where 

2 

q=u,d,s 


( 10 ) 


( 11 ) 


fS = |(> - E /")■ 


mNfn = {N\mgqq\N}. (12) 


2 
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2 
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lb 

2 

1 

V 

2b 

4 

2 + 2 

2 

4-h2 

S 

2b 

3 
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2b NLO 

- 

- 

5 

3 -h2 


2b NLO 

- 
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lb 

0-h2 

1 

2 + 2 

0 


lb 

2-h2 

1 +2 

A 

2b 

4-h2 

2 

2 + 2 

4 

P 

2b 

3 -h2 

5+2 


2b NLO 

- 

- 

5+2 

3 


2b NLO 

- 

A + 2 


Table 1: Left: leading chiral order of time (?) and space (x) components of the WIMP and nucleon cuiTents for vector and axial-vector interactions, for one-body 
(lb) and two-body (2b) operators. For the axial-vector nucleon operator, terms involving vertices from the NLO chiral Lagrangian (indicated by “2b NLO”) need 
to be included (see main text for details). The second number (“-I-2”) refers to the additional suppression originating from the NR expansion of the WIMP spinors, 
if momentum over WIMP mass is counted in the same way as for the nucleon mass. Right: leading chiral order of the WIMP and nucleon cuiTents for scalar and 
pseudoscalar interactions. 


Therefore, at leading order in the effect of integrating out the 
heavy quarks can be absorbed into a redefinition of Cg 




— y 


Q=c,h,t 


(13) 


For the u- and cZ-quarks the couplings are intimately related to 
the pion-nucleon cr-term <t„n I132r 


fN . 


rl\ _ 

Ju ~ 


2m\ 


A/;, /; = 


■N _ ^?rAf(l » fN 


2m\ 


(14) 


with ^ 


0.36 + 0.04 03311 and corrections A/" 


_ _ 

nifj+mii -i- yj.yj-r || -^-^| | l^j 

related to the strong proton-neutron mass difference via the 
low-energy constant cj. For the strange quark, the most ac¬ 
curate determination comes from lattice QCD ifj^l . The above 
Oja s) analy sis may not be accurate enough for the charm quark. 


see 023L135113611 for a study of higher orders in a^. 


This analysis generalizes to finite t if one defines 


mNfq(t) = {N(p')\mgqq\N(p)}, 
2 

= 


/g(0 


27 \ 


e^(t) = {N{p')\e>‘^\N(p)}, 


and replaces ^ f^(t) accordingly. 

The chiral expansion of cr„M starts with 


o-^N = -4ciM^ -h 0{py (16) 

in line with the O(p^) listed in Table [T] for the scalar one-body 
current. Note, however, that the power 2 does not imply a 
momentum-dependent coupling in this case, but a quark-mass 
suppression. As far as the f-dependence is concerned, the slope 
of the scalar form factors is dominated by nn scattering, which 
is known to not be adequately described by ChPT, but to require 
a reconstruction based on dispersion relations ll3T[^ . The t- 
dependence generated by other sources but light-quark scalar 
form factors was shown to be higher order in the chiral expan¬ 
sion in ifl^ . 


Defining 

Z /fw - 1 WG(0Cf), (17) 

q=u,d,s 

the NR one-body matrix element for the scalar channel be¬ 
come^ 

= xl'XrxlfN{t)Xs, (18) 

where (Xr'.s') are NR spinors for the incoming (outgoing) 
WIMP and nucleon, respectively. is of higher chiral or¬ 

der since the NR reduction of 75 produces a term -cr • q/(2m^), 
which we count as O(p^) for > niff. 

4.2. Vector 

The decomposition of the vector current at the quark level 
reads 

{N(p')\qrqmp)) = {N'lrFfit) - —q^FfitW), (19) 

2mff 

where the sign of the Pauli term is due to the convention in (HJl. 
To obtain a flavor decomposition of the vector current, one usu¬ 
ally assumes iso^in symmetry (corrections can again be calcu¬ 
lated in ChPT 1 I 41 I] ): 

F“'P{t) = Ff"(f), = F“’"(f), F'.'^it) = F/’"(f). 

' ( 20 ) 

In this way, one obtains 

F“’'’(f) = Ff’"(f) = 2Ff'^’'’(f) -H F™’"(f) H- F-yit), 

Ff’’{t) = F“’"(f) = F™’^(f) H- 2F™’"(f) H- F^it), (21) 

with electromagnetic form factors F™’^(t). At vanishing mo¬ 
mentum transfer this defines the vector couplings 

{Niqrqm = fyN\r\N), /; = - 2 /;^ = 2/^ = 2 . 

( 22 ) 


^The nucleon spinors include isospin indices according to xl^fNiOXs = 
5 ^ 5 ' [^/p(0 + /n(0)l + (fpif) ~ The Wilson coefficients match onto 

the conventions of O by means of the identification /az( 0) = "V^GfCo. 
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Corrections to (l22ll can be worked out in terms of magnetic mo¬ 
ments jiN - Qn + kn, electric radii as well as strangeness 

moments and radii {r\ , explicitly 



Ff^KN+0{t), Ff ^-KM-K^f,+0{t), 
Ff^K^^+0{t), 


corresponds to a second-class current i.e., it vio¬ 
lates G-parity, and will be ignored in the following. At vanish¬ 
ing momentum transfer only G^’^ contributes. Its coefficients 
are conventionally defined as 

{N(p)\qy‘y 5 q\N{p)} = Aq^iNly^jslN), (28) 

and isospin symmetry is assumed 

AuP = Ad", Au" = AdP, AsP = As". (29) 

The combinations 

- -a" - Au'’ - Ad'’ - gA, 

= Au^ + Ad^ - 2As^ ^3F-D, (30) 

are determined by the axial charge of the nucleon in the case of 
a 3 , or can be inferred from semileptonic hyperon decays for ag, 
yielding D xi 0.8, F « 0.46. The third combination 

(23) ^Au^ + Ad^ + As^ (31) 


with the Sachs form factors 


G^(f) = Ff (f) + ^F^(t) = gyv + 

4to^ 6 

Gj^(f) = Ff (f) + F^(t) = yuyvfl + + Oit^y (24) 


The NR one-body matrix elements involving a nucleon vec¬ 
tor current are 


^’^b.NR Xr'XrXs' 


/r(0 - ^ ■ (q - /tr X p)//’"'(f)^. 


+ ^ ‘ik'- ■ X q/2''’^(0T., 


2m^ 

^ . IT,, ,4 fA,Vy- 


“ 2m 


-Xff^Xr ■ 


—xy(T-’^Xrx'jr^t)xs 


2m, 


-T., 


(p - i(T X q)/f’'"(f) - i(T X q/^’'^(f)k,, 


(25) 


where 




q=u,d,s 


= ^ Z Cf'Ffit). (26) 


q=u,d,s 


4.3. Axial vector 


The decomposition of the axial-vector current at the quark 
level reads (see, e.g., Sil) 


<A(p')l^/r5?|A^(/^)> = (N'lrysGfit) - ys^Gf(t) 


2m\ 


'■^qvy5G‘!y\t)\N). 
2m M 


(27) 


is related to the spin structure function of the nucleon, it is not 
a scale-independent quantity. At = 5 GeV^ and G(q'j) the 
following values were obtained in 14511 

Au'’ = 0.842 + 0.012, Ad'’ = -0.427 + 0.013, 

As''=-0.085 + 0.018. (32) 

Besides the coefficients at zero also the momentum dependence 
of the flavor combinations 

0-4^ 1 _ 

= Qy^ys^Q = “ dyfjysd), 

1 

AI = Qy^ys^Q = — -pi^ytiysu + dy^ysd - 2sy^yss), (33) 
^ 2 V3 

can be analyzed in SU(Af) ChPT, but due to the anomalously 
broken U(1 )a current this is not the case for the isoscalar com¬ 
ponent. One obtains 

{N(p')\Almp)) = {N'\(ry5G\{t) - ys£^Gl{t)YA\N), 
{N{p')\Al\N{p)) = <A''|(/r5Gl(f) - r5^G^(f))i|A'), 


(34) 


with leading-order results 


Gi(t) = gA, Glit) 


3F-D 




G3(0.-^, gUo.-^J^. (35) 


f-M2- 


Empirically, the momentum dependence of G\{t), extracted 
from neutrino scattering off nucleons and charged-pion elec¬ 
troproduction, follows a dipole fit 


Glit) = 


{\-tiMir 


(36) 
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with mass parameter Ma around 1 GeV 142114311 . Since for gen¬ 
eral t the flavor structure cannot be inverted without additional 
input for the singlet component, we decompose the quark sum 
according tc0 

y CfGfpit) = Ct^G%{t) + C^^GlpiDr^ + Ct^G^pit), 


9 

with 


(37) 


“3 


V3 


'AA I /^AA 


cf ], 


^AA _ 

^3 ” 


^ r f^AA f^AA 1 

J’ 


cf = —[cf+Cf-2Cf], 


(38) 


and define 


^4,p(0 = ^[cf GXp( 0 + Cf + Cf Gl_/f)]. (39) 

In terms of these quantities, the NR amplitude reads 


Adf; 


,NR 


-xWXr -Xl 


o-8a(0 


4-m 


\o- ■ qg^(f) (40) 


Similarly, for the VA channel we define 


= ^[Co^G° ^(f) + C^'^Gl pit)? + Cl^G\ p{t)\, (41) 


to obtain 


= xlxrl^xl I O- ■ 0 - yy O- ■ q/*p(f) k 


^ 2m, 


In 


q P 


^friN 


(42) 




■y 


o-h’^(t) - —T^ ■ q^p(0 4rs- 


Am^ 

^rriN 


4.4. Pseudoscalar 

The pseudoscalar matrix element is usually parameterized as 
{N{p')\mqqiysq\N{p)) = {N'\mNG‘'f {t)ij=,\N). (43) 

By means of the Ward identity 

^ dpqy^y5q = ^ 2imqqysq - ^^Gf Gf, (44) 

9 9 

the corresponding form factor G‘^^’^{t) follows from G‘^^^(t) and 
G‘'p^{t), except for the singlet component, where the anomaly 
does not drop out. 


Gf f) = Gi,(f) + ^Gi,(f), 
4 ot^ 

Accordingly, we have 


/ = 3,8. 


= xl'Xr^xW ■ (t)xs, 

^^f,NR = ■ (iXr^xW ■ qh^{t)Xs, 


(45) 


(46) 


^At vanishing momentum transfer this equation maps onto the notation 
of by means of C^Aq^ = sflGfA^\{aQ + aiP). 


where 

85(0 = ^[crG^(0T' + cf G^ko], 

/iff) = A^C^PGl(t)T^ + Cg^^G^(f)]. (47) 


5. Two-body currents 


5.7. Scalar 

The scalar meson-exchange currents, involving both pion and 
q contributions, have been considered before in EES. The full 
expression reads 

\2 


M 


ss 

2 ,NR 


= fnMlxl'Xljl ■ T2 Xi2Xs,Xs2 

-xl'Xr[^) (48) 


where 


^ 0-1 ■ qi 0-2 ■ qz 

iql + Mf)iql + Mfy 


I = n, 7}, 


(49) 


,1 , ir.. 

pion decay constant F„ = 92.2 MeV j^ . Xsi (Xs') denote NR 
spinors for the incoming (outgoing) nucleons, with momenta p, 
ip'i)’ ‘Ji - P'i ~ Pi’ ^ - C/(77 -I- F), and the Wilson coefficients 
are collected in 


f ^ _L y C^S fn 

q=u,d 

with scalar meson couplings 
m,. 


fn 


J_ V f'i 
A 3 Zj « 


(50) 


f 

J U 


and 


■m 


= 0.32 + 0.03, /; = 


q=u,d,s 


m 

Mu + Md 


= 0.68 + 0.03, 
(51) 


r = - 

J ii o 


1 


Ml , 

4: = (6.9 + 0.4) X 10^^ 


r = - 

Jd n 


j__^ 

3 OT„ + Md 

md ^ (14 7 + 0.4) X 10^^ 


3 m„ + nid 


2 Mi. 

rq _ _ 

J8 - ^ 


■Ml.-Ml 


m2 

Lv±rj 


= 1.05. 


(52) 


One particular feature of the scalar two-body currents is that 
they cannot be written as a correction to the one-body coupling 
fy, since the scalar couplings of pions and n mesons probe a dif¬ 
ferent combination of Wilson coefficients OlOll . For this reason, 
even in the isospin limit they cannot be parameterized in terms 
of a single coupling cn as conventionally done for the one-body 
currents, see e.g. 1401] . 

5.2. Vector 

The only two-body vector current up to Oip^ appears in the 
AV channel 


M‘ 


,AV _ ^ f^Avl 8 a V 

' \2fJ 


2 ,NR “ A 2 2 

0-2 ■ q2 O"! 


ql + Ml 


xWxr-xlxli[Ti XT 2 Y 

^1 ^2 

+ (qi - q 2 )y 2 xs,xs2- 


0-1 ■ qi 0-2 


q? 


■Ml 

(53) 
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While the nucleon vector current itself has been studied in de- 
, the present application to direct detection is 


tail before 1118 
new. 


In fact, there are neither terms with i = 8 nor 77 contributions 
to i - 3. The reason for this can be traced back to the operator 
structure of the chiral Lagrangian: the coupling to the vector 
current occurs via a commutator of vector source and 

meson matrix. Expanded in Gell-Mann matrices, this leaves 
SU(3) structure factors and and the only non-trivial 
ones, apart from the direct couplings to the nucleon that led 
to drsl l. reduce to the SU(2) subset 


5.3. Axial vector 


6. Matching to NREFT 


NetO, we express our results in terms of the operator basis 


from 


il 


Oi = 1, 

O 4 = ■ s 


O 2 = (v■^)^ 


O 3 = iSn ■ (q X v-^). 


Ot -Sn y 
O[o - iSfi ■ q 


Os =/S;j, ■ (q X v-^), Oe =S^-qSivq, 
Og - iS^ ■ (Sn X q), 


Og = S;j- ■ 


Oil = ■ q, 

where S - <t/2 and the velocity is defined as 

K P 


2 my 2 m 


(56) 


(57) 


'N 


We find the relations 


The axial-vector two-body currents are 
1 


MiU = -X^X^ 




C6 

-trixq 




Fl 


2c,Ml- ** 


+ C3 


+ 2 d\T 


1^1 


q2 H- Ml 

/ _ 0-1 qq \ 

\ -<r Ml) 


0-2 ■ q2 


ql+Ml 


+ {l ^ 2 ) 

+ 2uf2[Tl X T2 ]^(o-i X tr2)(l - -<1 ^2 +Ml ) 


■XsiXs 2 , 


^l.NR ^ xWs'^^lfN^t^XrXs, 

^xlxlOiog3(t)xrX.s, 


^UNR = —Xl>xl 06 h‘^ {t)XrXs, 


o'2 ■ q2 

tn^ 


\ql + Ml 




niN 


1 




(f04 H- Ofl)fl[’^{t'^rXs 


^f,NR = Xlxl 208/f•'^(f) -H —^9(/f’^(0 + f 2 ^{t)) XrXs, 


-Ar'As' 

xtr'A^' 


- 404g"(f) -I- -^ 0 (,g^p(t) XrXs, 


— '101 + — 0c^ 


hy.i)XrXs- 


(58) 


(54) 


where the terms that do not contain an explicit q-dependence 
(q = -qi - q 2 ) and the ce-term are taken from MM, while the 
finite-q pion-pole corrections were derived in ||l5|] . The AA 
two-body current as in 02111 has been applied in the calculation 
of structure factors for spin-dependent scattering in lilllil, 
whereas the two-body current in the VA channel. 


1 _ t t . .t I .•!_ _ n3 O'! ■ q2 


A^I^nr = --2Cy^2y>XrXl,xl\i[r, Xr2] , 


^2 3 2/^2xvr'xkA/t,y'/v 

+ (1 ^ 2 )\x^iXs 2 , 


ql + Ml 


(55) 


has not been considered before. 

For similar reasons as in the vector case there are no / = 8 
or 77 contributions from the leading-order Lagrangian. In prin¬ 
ciple, one could calculate corrections from the NLO SU(3) 
Lagrangian, in analogy to the SU(2) result for AIjnr- How¬ 
ever, there is a lar ge n umber of poorly-known low-energy con¬ 
stants (see 0 or l48l] for the matching to SU(2)), which would 
severely limit the predictive power. 

Finally, due to the derivative in the Ward identity (l44l l. there 
are no pseudoscalar two-body currents at Oip^). 


This shows that as a result of QCD effects, the operators in 
the NREFT are not independent. For example, both axial and 
pseudoscalar operators combine in the nuclear matrix element 
In addition, up to 0{p^) only 8 of the 11 operators 
of (l56l l are present. However, because itself enters only 

at 0{p‘^), they are mapped onto 7 amplitudes, so that the re¬ 
lations cannot be inverted. This is because and 

involve the three operators Og-g. This implies that some opera¬ 
tors, e.g. d?6, can be isolated by having a particular quark-level 
interaction, but this is not possible in general, as demonstrated 
by the example of Og-g. If we retain subleading corrections 
in the NR expansion of the spinors, the missing operators ap¬ 
pear, accompanied by additional combinations: O,, in terms 
of Af™ R’ <^2 and 0 $ in O 3 O 3 in and OiO% in 

A 4 AA 

2’'‘i,NR' 

In the limit where becomes (significantly) larger than the 
nucleon mass also should be dropped, as well as the 

1 / 777 ^ suppressed terms in and Adj^^j^. In contrast, all 

two-body currents up to 0{p^) are independent of 777 ^. They 
appear in the SS,AV, AA, and VA channels. 

We stress that the above discussion merely pertains to the 
mapping of operator structures, it does not take into account the 
evolution of the scale dependence that is required when match¬ 
ing the coefficients of a pionless theory, valid for scales below 


6 




























the pion mass, and ChEFT, defined at chiral scales. This in¬ 
volves also effects related to the limitations of the “Weinberg” 
counting scheme applied here ifi^ . and would have to be taken 
into account in the matching relations required for translating 
NREFT coefficients to the QCD scale. In addition, there may 
be effects from operator mixing, originating from the interplay 
between the nucleon-spin dependence in the ChEFT WIMP- 
nucleon scattering operator and that in the high-momentum part 
of the ChEFT NN potential, which would also have to be con¬ 
sidered when evolving NREFT operators to the QCD scale. 

7. Summary and discussion 

In this Letter, we have developed the constraints that chiral 
symmetry of QCD imposes on the nuclear matrix elements that 
can enter in dark matter direct detection. We provide explicit 
expressions for one- and two-body currents in WIMP-nucleus 
scattering for vector, axial-vector, scalar, and pseudoscalar in¬ 
teractions up to third order in the chiral expansion. The chiral 
power counting, summarized in Table[T] shows that at this order 
there are two-body currents that have not been considered and 
may be of similar or greater importance than some of the one- 
body operators, see (l5^ and (l55l l. Moreover, the matching to 
NREFT shows that not all allowed one-body operators appear 
at this chiral order and that the operators in the NREFT are not 
independent. 

The chiral power counting applies to the one- and two- 
nucleon level. In nuclei, the different interactions can lead to 
a coherent response that scales with the number of nucleons in 
the nucleus or to a single-particle-like response. In a next step, 
we will evaluate the nuclear structure factors, including the con¬ 
tributions from two-body currents, and provide a set of response 
functions for the analysis of direct-detection experiments. This 
will also allow us to assess how constructive or destructive the 
interference of operators based on the constraints provided by 
chiral symmetry proves to be. 
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